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^ ' Abstract. We studied theoretically a coherent population trapping resonance 

O^. formation in cylindrical cell without buffer gas irradiated by a narrow laser beam. We 

take into account non-zero probabilities of elastic ( "specular" ) and inelastic ( "sticking" ) 
^S] ■ collision between the atom and the cell wall. We have developed a theoretical model 

^ I based on averaging over the random Ramsey pulse sequences of times that atom spent 

in and out of the beam. It is shown that the shape of coherent population trapping 
resonance line depends on the probability of elastic collision. 
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1. Introduction 

In the simplest case, the coherent population trapping (CPT) effect appears in three- 
level A-system, see figure [1] The atom is interacting with a two-frequency laser field 
coupling two metastable states with a short-living state. When the frequency difference 
between the laser field components is equal to the frequency of the transition between 
two metastable states (two-photon resonance conditions), the "dark state" that does not 
interact with the laser field appears. This "dark state" is the coherent superposition 
of the metastable states. If one scans the frequency difference of two laser components 
near the two-photon resonance, one observes the CPT resonance, an abrupt decrease 
of the fluorescence intensity [lH3]. CPT phenomenon allows to create a window in 
the absorption spectrum [1]. This effect is used in compact frequency standard OE], 
the stimulated adiabatic passage technique of population transfer between the states of 
atoms and molecules [7], laser cooling [8J and magnetic fleld measurements J9|[T0]. 

It is well known [TTHT9] that the CPT resonance shape is affected by the movement 
of atoms between the beam zone illuminated by the laser radiation and the dark zone. 
The theory of the CPT resonance formation in wall-coated cells without buffer gas based 
on averaging over random Ramsey sequences of times which the active atom spent in 
the bright and dark zones was developed in [T214H] . In these models, it was assumed 
that in any atomic-wall collision the atom sticks to the cell wall for some time. During 
this time atom exchanges its kinetic energy with the cell wall and then it is released back 
to the cell volume with alternated velocity that is not correlated with atomic velocity 
before the collision. Atoms make many transitions between these two zones during the 
lifetime of metastable states coherence. This cohrerence is assumed to be preserved in 
the collision [20] . 

Assumption about sticking collisions is based on the studies of relaxation of 
polarization in Rubidium vapors in paraffin-coated cell performed by Bouchiat and 
Brossel [21J. These studies showed that at least some atoms sticks to the cell wall for 
some time in the wall collision act. During this time, atom exchanges its kinetic energy 
with the cell wall and then it is released back to the cell volume with alternated velocity. 
At the same time, some part of atoms can collide elastically with the cell wall. 

In this work, we make the generalization of the model [T2Hllj to the case when there 
is a nonzero probability a (refiection coefficient) of elastic collision between the atom 
and the cell wall. Examples of an elastic collision are evanescent wave mirrors [22|[23] or 
magnetic gradient mirrors [23H26] for cold atoms. We assume that an elastic collision 
do not change velocity components parallel to the cell wall (normal component changes 
sign). We study the infiuence of the refiection coefficient a on the CPT resonance 
line shape at different intensities of laser radiation provided that the width of the laser 
radiation is much smaller than the Doppler width of the optical transition. 

The density matrix equations describing the atomic polarization inside and outside 
the laser beam are presented in section O In section [HI we discuss the influence of elastic 
and inelastic collisions on the atomic motion in the cell and derive the expression for 
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the excited state population. The procedure of the ensemble averaging is described in 
section HJ Results of calculations are presented in section [51 Conclusions are made in 
section |6l 

2. Basic equations 

Let us consider a gas of three- level atoms with the excited state |3) and metastable lower 
states |1) and |2) (one of which can be the ground state). The field component with 
frequency uji couples the states |1) and |3), whereas the component with frequency u^ 
couples the states |2) and |3) (A-scheme of the atom- field interaction, see figured]). The 
interaction of the field with the atoms is described by Rabi frequencies Vi = |(ii3 ■ Ei/h, 
V2 = ^c^23 ■ Ei/h, respectively. Here di3 and ^23 are the matrix elements of the atom's 
dipole momentum, Ei and E2 are the amplitudes of the laser field components. The 
one-photon frequency detuning of the first component is Ql = uji — oou, whereas 
Q = Ui — U2 — UJ12 is the two-photon detuning. Here, Uij is the transition frequency 
between i and j states. In our case 1^12 belongs to the microwave range. 
The equation for density matrix describing the atom in the field reads 

Pij — ~T /_^ Wik Pk'j — PikHkj] + 2_^ '^ijM Pkl ; (1) 

k k,l 

where Hik are the matrix elements of the Hamiltonian H = Hq + hV{vz,t), Hq is the 
Hamiltonian of the free atom. The operator of dipole interaction of the atom with the 
laser field HV{vz, t) depends on time t and on the component of atomic velocity Vz along 
the field propagation direction; Fj^ ^^ are the elements of the relaxation matrix. The 
nonzero matrix elements are 

ri2,i2 = r2i,2i = — r^ 

rii,22 = r22,ii = — rii,ii = —1^22,22 = r, 

TlS.lS = r31,31 = r23,23 = 1^32,32 = "7 ) (2) 

ril,33 = r22,33 = 7/2, 
r33,33 = - 7- 

Here 7' = (7 + ^l)/'^ is the relaxation rate for optical coherences pis, P23 ^27j, Tl is 
the spectral width of the laser radiation, 7 is the spontaneous relaxation rate for the 
excited state, F is the relaxation rate in the ground state. 

The light absorption in the cell is proportional to the excited state population P33. 
We consider the case of weak fields V^i.2 ^ 7- In this case P33 is much smaller than 
the populations pn and P22 and can be expressed via pn, P22 and pi2 using standard 
adiabatic elimination procedure [28] . 

Usually this procedure is applied in the case of large detuning Q^ ^ V^i 2, but it 
is also valid in the case of 7 S> Vi,2- Let us consider for illustration an open A-system, 
when the atoms due to spontaneous emission decay from the excited state |3) to the 
states, different from |1) and |2). It is well known that the irreversible loss can be 
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Figure 1. A-scheme of the interaction of ttie three- level atom with the field, Vi and V2 
are the Rabi frequencies, T and 7 are relaxation rates for ground and excited states, 
respectively, and Vt^ and 17 are the one-photon and two photon (Raman) detunings. 



incorporated in the Hamiltonian [29] . For the open A-system in the rotating frame the 
Hamiltonian H is 

/ Vi \ 

H = h -n V2 . (3) 

\ Vi V2 -VIl - 11/2 j 
The set of equations for the amphtudes is 

(4) 
C3. 

Typical evolution time of Ci and C2 is about 1/V^i,2- If I^l — i7/2| ^ Vi^25 "we can 
put Cs = in dl]) on the time scales of evolution of Ci, C2 and then express C3 via 
Ci and C2, i.e., to perform the adiabatic elimination of the excited state. In the case 
of closed A system considered in our paper, the Hamiltonian H is the Hermitian part 
of H. Considerations for the density matrix equations similar to ones from previous 
paragraph can be performed analogically. Another important point is that in the case 
considered in this paper (room temperature and laser beam diameter of order of 1 mm 
or higher), the typical time of flight of the atom through the beam or dark zone is much 
larger than I/7'. 

We also neglect the Doppler shift of the frequency of microwave transition or, what 
is the same, we suppose that the Doppler shift is equal for two optical transitions. 
This approximation is vahd when the relative phase between two optical fields doesn't 
changes significantly on the length of typical path of the atom (Dicke narrowing [30]). 
In [3HI32]. it was shown that for the longitudinal cell size smaller than A21/4 (where 
A21 = 27rc/ci;2i is the wavelength of transition between the states |1) and |2)), the shape 
of CPT resonance in coated cell without buffer gas practically does not depend on the 
cell length. Therefore we can neglect the Doppler broadening of the microwave transition 
if the cell size is less than A21/4 that is equal to 1.1 cm for microwave transition in ^^Rb 
atom. 
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Using the normalization condition 

Pii + P22 = 1 (5) 

and denoting p = {f,R, J}, where f = pu — P22, R = R-e(pi2), J = Ini(pi2), we can 
rewrite the equations for the density matrix for the ensemble of atoms in the laser beam 
in the form: 

/ = G^^^V^ - (ly + r)/ - 4F^ J , 

R= - G^ -iW + T)R-{n-A)J, (6) 

7' 

j = F^^ f + (n - A)R- (W + T)J . 

7 

Here G = G{vz) and F = F{vz) denote the real and the imaginary parts of the expression 
1' /il' ~i{^L — kvz)) respectively, k = ui/c is the wave number of the optical radiation, 
A = F(V/ - Vi)/^ is the light shift, W = GiV^ + ¥2)/^ denotes the optical pumping 
rate. The excited state population expressed in terms of / and R is 



W G 

P33 = — + — J 



{y^-V^)-f + W^^2-R\. (7) 

To simplify the calculation procedure we use the matrix notation formalism similar 
to [35. The set of equations (|6]) can be symbolically presented as 

'p{Vz,t)=k{Vz)p{Vz,t)+B{Vz), (8) 

Also we can write the expression (ITj) in the form 

P33(t^., t) = U^{Vz)p{Vz, t) + V{Vz) . (9) 

Furthermore, for sake of brevity we will skip the argument Vz in A, B, U and V . 

The evolution equation for the density matrix of atoms outside the laser beam can 
be obtained from equations ([6]) by substitution Vi = V2 = W = A. = Q. One can rewrite 
the resulting set of equations symbolically as 

'p{Vz,t)=Mp{Vz,t). (10) 

Note that A' does not depend on Vz- 

Finally, the density matrix of the atom is described by the equation ([8]) inside the 
laser beam and by the equation flTOj) outside the beam. The solution of the equation ([8]) 
can be presented as 

p{t) =(l- e^-(*-*«))p5 + e^-(*-*«)/)(to) , (11) 



where ps = —A^^B is the stationary solution of ([8]) and to is the initial time for which 
the density matrix is known, I is the identity matrix. The solution of equation (ITOl) is 
written in a similar way: 

m = e^'-(*-*")p(to) . (12) 
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Figure 2. The trajectory of a single atom in a cell in the beam passing regime (a) and 
in the dark regime (b). (c) dependence of the optical pumping rate W on time. The 
magnitude of W depends on the velocity of the atom that causes the Doppler shift of 
the atom's transition frequencies. The checkmarks show the times of elastic atom-wall 
collisions and the asterisks show the times of inelastic collisions. 



3. Atomic motion in the cell 

In this paper, we consider a cylindric cell illuminated by a co-axial cylindric laser beam. 
Let us imagine first that the atom undergoes only elastic collisions with the cell walls. 
Therefore there are only two possibilities: either the atom passes through the beam 
zone during each pass through the cell as it is shown in figure [2]^a), or it does not enter 
the beam zone at all (figure [2]^b)). In the first case we say that the atom is in the beam 
passing regime^ while in the second case the atom is in the dark regime. Obviously, the 
resonance is formed only by the atoms in the beam passing regime. Let us consider some 
of such atoms. At the observation time tA this atom was in the beam zone during some 
time t continiously. Before the enter to the beam zone the atom was in the dark zone 
during time r' . Before this the atom was in the beam during time r. In turn, before 
this the atom was in the dark zone during time r' etc. Using equation (ITT]) and ( !T2|l . 
we obtain 



pit 



A 



I-e" 
+ e^-^e'" ■ 
I-e^-* 



Ps + e 

I- 



+ e" 



I- 
Ps + e"* 

I-e^' 



Ps 



I 



Ps- 



(13) 



Now let us consider more common case when some of the wall collisions change 
atomic velocity. After such a collision the atom being initially in the dark regime can 
either stay in this regime (which does not affect the equations for the density matrix) 
or switch to the beam passing regime. An atom that was being in the beam passing 
regime before the collision can either pass into the dark regime or get back to the beam 
passing regime, but with different value of the velocity and angle of reflection with the 
wall of the cell (and therefore with different values of r and r'). For sake of definiteness, 
let us call the time of entry of the atom into the beam zone after an inelastic collision 
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as the beginning of the beam passing regime, and the time of exit from the beam zone 
as the end, respectively. Obviously, if the atom undergoes N elastic collisions with the 
cell wall in the beam passing regime, then it goes A^ times through the dark zone and 
A^ + 1 times through the beam zone (see figure [2]^c)). 

The density matrix p^ describing the atom at the time of exit from the beam passing 
regime can be obtained similarly to the expression flT3|) . It reads 



Pe 



I + e^-^e^ ■"+...+ ( e^-"e^-^ 

N 



N- 



{l-^-^)ps 



+ (e^-"e^-" ) e^-"pb 

N 



I - e^-^e 



N+l 



(e^-e^'-') ](l-e^-)p5 



+ e-^-'e^ 



e^-^Pb. (14) 

Here, pb is the density matrix of the atom at the time of beginning of the beam passing 
regime. It can be expressed via the density matrix Pe,-i for the previous time of exit from 
the beam passing regime (it is indicated by the index "—1") using the expression ( IT2l) 

Pb = e'-'-^'Pe,-! , (15) 

where r^ is the time that atom spent in the dark. Similar to (fT^ . one can obtain the 
expression for the density matrix at the observation time tA'- 



I 



p(t^)=e^-*e^-^(^I-e^-^e^-^ 

+ (l-e^-*)p5 + e^-*(e^'-^V--)"p,. 
The expression for the excited state population P33 reads 



(^^K■r^K'■ry^^(^_^A■ry^ 



(16) 



tT 



P33 = f/Me"-*e^-^ I-e 



I- e^-V 



l-^-^]ps 



+ (I-e^-Mp5 + e^-*(e^'-^'e- 



Pb\ + V. 



(17) 



Here, n is the total number of collisions in the beam passing regime immediately 
prior to observation time Ia- Averaging the expressions (fT7|l . (fT^ and ( !T5|) over all 
the atoms in the beam zone, we obtain the average excited state population (pss) that 
determines the absorption of light in the cell. 



4. Ensemble averaging 

4.1. General scheme of averaging 

Let us consider in detail the averaging of f lT7|) . In this paper, we consider the excitation 
of the CPT-resonance by a laser whose spectral width is much smaller than the Doppler 
width of the optical transition. Hence, the main contribution to the resonance is given 
by the atoms with small longitudinal velocity, that are rarely colliding with the ends of 
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the cell. We assume that the collisions of the atoms in the beam passing regime with 
the side walls of the cell do not change Vz- The average population (P33) is 



(P33) 



+ {V)+ ''" 
X I — ( e" ■ e 



U'^i 



(p.) 



U^e' 



I 



Ps 



(18) 



where the angle brackets denote the averaging over the z-component of the velocity v^, 
over the times t, r, r', Td and over the numbers n and A^ of elastic collisions in the beam 
passing regime prior the observation time Ia and prior the exit from the beam passing 
regime (in the expression for {pb)) respectively. The averaged density matrix (pt) of the 
atom entering the beam passing regime can be obtained by averaging of expressions ^\M 
and (IIS]): 
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C 



l-e' 



-1 



7V+1- 



I-e^ 



Ps 



N 



(19) 
(20) 

(21) 

(22) 
(23) 



Let us average over n and A^ first. Assuming that the probability of the elastic 
scattering a does not depend on the velocity of an atom and its internal state, we 
obtain the probability P{n) = (1 — a)a^ for n consecutive elastic collisions. Hence, for 
any matrix R we get the matrix R" averaged over n as 



(R")„ = X]P( 



n] 



;i-a)(I-aM) 



-1 



(24) 



n=0 



Exactly the same kind of distribution P{N) = (1 — a)a^ is valid for a total number of 
elastic collisions A^ in the beam passing regime; therefore, the matrix R^ averaged over 
A^ is described by the same expression: (R^)jy = (R")„. 

The next steps are the averaging over times t, r, r' and t^. It is obvious that the 
averaging over the time Td of atomic stay in the dark regime is independent on t, r, 
r', while r and r', in general, are determined by transversal component of the atomic 
velocity in the beam passing regime and by the angle between the velocity and the 
wall of the cell. The time t is determined also by atom's position in the beam zone at 
the observation time t^. Strictly speaking, one must integrate over four variables t, r, 
r' and Vz- However, we make a simplifying assumption which reduces the number of 
integration variables to 1. 
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In the beam passing regime atom goes A^ times through the dark zone and A^ + 1 
times through the beam zone. Each passage through the dark zone takes time r' and 
each passage through the beam zone takes time r. Obviously, for each passage the times 
t' and T are not independent. According to our assumption we consider times of different 
passages of the dark and hght zones as independent random variables. We emphasize 
that this assumption refers only to ensemble averaging. In this case the averaged values 
of e*"^, e*'"^, e* "^ and e^ "^"^ in flT8|) - fl23|) can be calculated independently. Consider, 
for example, the averaged expression of e"^'^. Its easy to present this expression in the 
form li] 



(e^-), = X 



/ (e^^^)^ 

(e^^Or I ^'\ (25) 

\ (e^3^), 

where X = {X^^\ X'^'^\ X^^^) is the matrix constructed of the eigenvectors of A 
corresponding to the eigenvalues Ai,A2,A3. Therefore, averaging over r is reduced to 
the calculation of integrals 



(e^Or = <7.(A) = / e^--^(r)dr, (26) 

that can be performed analytically. Here, fr is the probability density function of the 
random variable r. Averaging over t, r' and r^ is analogous to (1261) with functions 5't'(A), 
gt{\) and gr^^X). Note that Re(A) < 0. 

4-2. Averaging over t and t 

First of all, we find the probability density functions f^ and ft of the random variables r 
and t. The time r that the atom spend in the beam zone is 



2^r^-Wsin\ 
T = , (27) 

v± 

where ip is the angle between the atomic velocity projection on the plane orthogonal 
to the cell axis and the normal vector to the cell surface, R is the radius of the 
cell, r is the radius of the beam and v± is the component of the velocity orthogonal 
to z. Probability density function of the transversal velocity v± is the two-dimensional 
Maxwell distribution M2{v±) = 2exp{—v'j_/v^)v±/v^ where vt = \/2kBT/m is the most 
probable atomic velocity, ks is the Boltzmann constant, T is the temperature, m is the 
mass of the atom. The probability density $(</?) of the random variable ip is proportional 
to cosy? |U1[M1[35]. Taking into account that in the beam passing regime the p takes 
values between and arcsin(r/i?), we obtain $(</?) = Rcosp/r. Then, the cumulative 
distribution function -Ft(t) of the random variable r is 

arcsin(r/i?) oo 

Fr{r)= U{v) I M^Mdv^dp 

2y' r'^-R^ sin^ ip 
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and the probability density function is 

dFJr) 



fr{r) 



1 
- + 

r 



4r 

e " "t 



W---T^,rf 



2ir 



2„2 



r^f. 



dr r 4ir Vt^t 

Note that the average time that the atom moves through the beam zone is 

7^3/2 ^ 



/•OO /'OO 

/ /(r)rdr= / (1 - F(r)) dr 



2 t;T 



(28) 



(29) 



(30) 



Let us find the probabihty distribution function ft of time t that atom stays 

continuously in the beam zone before the observation time Ia- The averaging is carried 

out over all the atoms in the beam zone. Let dn/dt atoms arrive to the beam zone per 

unit of time. Then the total number of atoms in the beam zone is 

dn r°^, _ 

(31) 



N, 



tot 



dt 



POO 

/ {I - F,{t')) dt' . 
Jo 



Here, we took into account that 1 — F^{t') is the atom's probability to remain in the 
beam zone after time t' from the entrance to the beam. Therefore the number of atoms 
that have already spent time t or more in the beam zone is 






•*+ 



F,{t'))dt'. 



Hence, the cumulative distribution function Ft{t) of the random variable t reads 

j^{l-F^{f))dt' 



m) 



r^i 



;o v^ FAt'))dt'' 
and the corresponding probability density function 

1 _ F,{t) 1 - F,{t) 



IS 



Mt) 



m 



/;°(i-F.(to)dt' 



T 



TVt^TX -^i^ 



T erf 



2ir 



(32) 



(33) 



(34) 



(35) 



4ir \ TVt , 

Using (134|) it is easy to express 5'r(A) in terms of (?t(A): 

gr{X) = \^\fgt{X) (36) 

Note that the function 5't(A) = J^ e'^^ ft(t) dt cannot be calculated in terms of 
elementary functions. We can find approximate but close to accurate result by replacing 
the precise expression 03 5 p by some approximation similarly to 

/,(*) -* ^/, (*^^ 

r V r / 
2 



Mx) 



-bx 16(1 
:e + 



,— az\3 



+ Axe 



(37) 



7]-3/2 3y^3/22;2 

We choose the values of coefficients A = 0.70839, a = 0.49, b = 2.286 and c = 1.272. 
These coefficients ensure (1) the same asymptotical behaviour of ft{t) and ft(tvT/r)vT/r 
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Figure 3. Plots of ft {x) (dashed line) and r-ft{x-r/vT)/vT (solid line) 



at infinity, (2) coincidence of ft{t) and ft{tvT/'r)vT/T near the zero to the second order 
with respect to t and (3) correct normahzation j^ ft{x) dx = 1. The plots of the 
functions r ■ ft{x ■ r /vt)/vt and /t(a;) are presented in figure [3l 

The approximation fl37|) makes it possible to find the analytical expression for gt{X): 

9tW = gtiXr/vr) , 



~9M) 



2 1 A 

+ 



7T^/H-X (c-A)2 

16 

+ ^^^[P(-^) - 3p(a - A) + Sp{2a - A) - p(3a - A)] , 



where we denoted 

p{f3) = f3\nf3. 
It is easy to find (7r(A) from fl38l) and fl36|) . 



(3^ 



(39) 



4-3. Averaging over t' 

Now let us calculate the function 5't'(A). The time r' that the atom being in the beam 
passing regime spend in the dark zone continuously is 

, 2 



r = — \ R cos (y9 — y r^ — i?^ sin (f 



(40) 



If the beam diameter is much smaller than the cell diameter, the expression in brackets 
can be approximately replaces by the constant 

r"arcsin(r/_R) 



R COS (p — y r"^ — R'^ sin (f \^((p)d(p 



r^ vrr R r 

i?A/l-^-y + -arcsin- 



(41) 
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that equals to average transversal path length in the dark zone. This approximation 
allows to readily obtain the cumulative distribution function F' of the random 
variable r': 



F'ir') 



v^ 



and 



Vx 



(42) 



(43) 



An attempt to directly calculate the function 5't-'(A) = X f^{l — F' {x))e^^ dx does not 
yield the result that can be expressed through elementary functions. As in previous 
case, we can find approximate but close to accurate result by replacing the precise 
expression (H2]) by 



F'(r') -^ F' ( r'^ 



F'(x) 



(l_e-<^(^-^o)) 



-^(^-^o)(l + 6(x-Xo)), X>Xo 

X < Xq 



(44) 



The values of the parameters a = 0.89, b = 2.56 and xq = 0.3864 ensure smoothness 
of the function F'{x) anywhere, good approximation of F'{xi±/vT) in the region where 
it is appreciably different from zero (see figure H]) and get the correct value of time f'. 
Using the approximation function ( l44l) . we find 

gr'{\)^~g'{\i±/vT), 



g'{A) = Ae^^'o 



26 -A 
[{b^Tij 



2 + p(4a - A) - 4p{3a - A) 



-4p{a - A) + 6p{2a - A) + p{-A) + - 



(45) 




Figure 4. Plots of functions F'{x) (daslied line) and F'{x£±/vt) (solid line). 
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4.4- Averaging over Tci 

The last function we should find for calculating the fiuorescence of the atoms in the cell 
is Qr^^X). The exact calculation of the probability density function fd of the random 
variable r^ is a nontrivial problem, because in the dark regime the trajectory of the 
atom consists of a random number of segments of random length; along each of them 
the atom moves with a random velocity. Nevertheless, it is clear that the probability 
density exponentially tends to zero when t^ — ?■ 0. The minimal distance that an atom 
can pass in dark regime in the transverse direction is 2{R—r). The transverse velocity v± 
of the atom with the probability greater than 98% is less than 2vt, so we assume that 
fd{Td) = when Td < tq = {R — r)/vT- If the beam diameter is much smaller than the 
diameter of the cell, the atom usually returns from the dark regime into the beam passing 
regime after many wall collisions. The probability for atom to return to the beam passing 
regime after certain collision does not depend on the number of previous collisions. Such 
processes are usually described by an exponential probability distribution function, so 
for the function fdi^d) of the random variable r^ we take 

fd{rd) =< (46) 

[ U , Td < Tq, 

where h = l/(frf — tq), r^ is the average time of atom's being in the dark zone. From (126|1 
and (HED it is easy to get gr^{X) = e^''°/{l - X/h). 

The value of f^ can be found from the "time balance" conditions. It means that 
the total time that an atom spends in the beam or in the dark zone is proportional 
to the volume of this zone. Consider the cycle in which an atom being in the beam 
passing regime undergoes N elastic collisions with the wall, and then moves into the 
dark regime. The history of any atom can be represented as a sequence of such cycles. 
The total time tf, spent by an atom in the beam zone over the cycle is t^ = {N + l)r. 
In turn, the time td spent by the atom in a dark zone during the cycle is td = Td + Nt'. 
The ratio of averaged times tf, and td over the different cycles must be equal to the ratio 
of volumes of the bright and dark zones. Hence, we find 

where N = a/(l — a) is the average number of successive elastic collisions with the cell 
wall. 

In the conclusion of this section we write down the expression for the averaged 

excited state population: 

00 

(P33)= [ M^{v,)U^Gt{l-a)(l-aG'^,G'^,y'dv,- (p,) 



—00 
00 



+ / M,{v,)Vdv,+ / M,{v,)U 



(I - Gt) + GiC;, (4^ 
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X (I - g,g;,)"^ [i - (1 - tt) (i - «g,g;,) 



Psdv^ 



G, =X 



G^ =X 



where Miivz) = e ''""'^^/(i'tV^) is one-dimensional Maxwell distribution function, 

/ ^t(Ai) \ 

gt{\2) X-\ 

\ gtih) J 

( 9r{\l) \ 

gr{\2) X-\ 

V g^{\^) I 

I qAK) \ 



g:, =x' 



gr>{\'2. 







V 











X' 



/-I 



9r'{>^z) ) 



(p,) = (^i - (g;;dj (g;^)c, 

where expressions f l22|23l) for C and D can be rewritten as 



(49) 

(50) 

(51) 
(52) 



C 



M,{v,)(l 



X (l-aG^G;, 



-1 



I-(l-a) 



(I -Qr)Ps dVz 



D 



/ Mi(t;^) G^(l -a){l- aG;,G^) dt;^ 



') 



Td 



M,iv,)X' 



( 9rAK) \ 

g^X',) 

V ^-.(^3) / 



X'-^ dv. 



(53) 
(54) 

(55) 



5. Results of calculations. 



An absorption of the laser power in the cell filled by the gas of three-level A-atoms 
is proportional to the average excited state population described by expression (fTSl) . 
Typical parameters of ^^Rb atom in vacuum cell were taken as the parameters of the 
A-atom in our calculations. The atomic mass m was assumed to be equal to the mass of 
^■^Rb isotope, the temperature T = 20°C, the ground state relaxation rate T = 300 s~^, 
the optical coherence relaxation rate 7' = 1.8 x 10^ s^^. We assume V^i = V2 = V^ 
and we set the optical detuning Ql = 0. Calculations were performed for the cell 
radius R = 0.5 cm and for three values of the laser beam radius r {r = 0.5, 1.5 and 
5 mm) and for five values of the probability a of the atomic-wall elastic collision {a = 0, 
0.25, 0.5, 0.75 and 1). 

The structure of the CPT resonance for a = was discussed in pTHT^ . Resonance 
consists of a broad pedestal whose width is about several tens of kilohertz and of a narrow 
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Figure 5. The shape of CPT resonance for different a, V and r. Parameters: 
r = 1.5 mm, F = 7.7 x 10^ s^^ (a, b); r = 0.5 mm, V = 1.46 x lO'^s"^ (c); r = 0.5 mm, 
F = 7.3xl05s-i (d); r = 5mm, F = 3.3x10^ s^^ (e); r = 5mm, T^ = 1.46xl0'^s-i (f); 



central peak (see figure l^a)). Broad background is formed by tlie atoms tliat liave 
passed tlirougli tfie beam zone only once, whereas the central peak is formed by multiple 
passes. If the elastic collision probability a is nonzero, the additional "intermediate" 
structure in CPT resonance appears. This is a peak whose width is greater than the 
width of the narrow central peak, but less than the width of the broad background. This 
peak is formed by atoms that pass the beam zone many times in beam passing regime. 
It should be noted that in the elastic collision act the atomic longitudinal velocity Vz 
remains constant, in contrast with the case of inelastic scattering when atoms pass the 
beam zone just once in the beam passing regime. This fact explains the occurrence 
of "intermediate" structure depending on a when the cell radius coincides with the 
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radius of the beam (see figure O^e) and (f)). Indeed, in the considered case of the laser 
field with narrow spectral linewidth only the atoms with small |f^| contribute to the 
CPT resonance. The "intermediate" structure is formed by the atoms that remains in 
this velocity group during several passage through the cell volume whereas the narrow 
structure is formed by the atoms that change their velocities after one wall collision. 

At Fig. 5(a)-(f ) we see that the higher elastic collision probability a corresponds to 
a lower excited state population at fi = 0. To explain it we should note that for higher 
a atoms spent longer time in the beam passing regime; therefore, they have more time 
to came into the dark state. 

It should be noted that if a = 1, the narrow peak disappears. It can be explained 
by the fact that all the atoms contributing to the CPT resonance formation are in the 
beam passing regime permanently. These atoms spent relative higher part of time in the 
beam zone than all the atoms in case of a = 0. Therefore, the narrow peak is formed by 
atoms which passed the beam, then spend quite a long time in the dark regime and then 
returned back to the beam. Intermediate structure is formed by atoms which passed 
the beam several times repeatedly in the beam passing regime. The time between two 
subsequent beam passage occurs to be smaller than in the case of inelastic collision. 
Therefore, the resonance occurs to be broader. The broad pedestal is formed by atoms 
that pass the beam zone once during the coherence time. 

6. Conclusions 

We developed a theory of the coherent population trapping resonance in a cylindrical 
cell filled with antirelaxation coated walls. The laser beam is assumed to have constant 
intensity in cylindrical region. The atom-wall interaction is described by the probability 
of elastic collisions. It is shown that the nonzero probability of elastic collision leads 
to distortion of the CPT resonance lines that consists in appearance of the additional 
peak formed by the atoms, repeatedly passing through the beam zone in the beam 
passing regime. The general scheme of calculation can be applied not only to CPT 
formation in wall-coated cells but also for the cell with small buffer gas pressure (much 
lower than the necessary for the diffusion approximation applicability). Another possible 
application can be connected with the high-precision clock based on an ensemble of cold 
atoms or ions, e.g. ions of ^^^Th [37]. The energy splitting of the ^^^Th ground-state 
doublet is about 7.6eV [38] corresponding to vacuum ultraviolet and long, about an 
hour, expected half-life of the excited state of the nuclei [39] makes this nuclei a possible 
promising object for nuclear frequency standard. The study of one-photon resonances 
in the two-dimensional trap with zonal pumping may be applicable to the spectroscopy 
of this nucleus. 
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